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1. Introduction 

Heavy-quark photoproduction is a phenomenon of considerable interest. It is 
closely related to the hadroproduction phenomenon, but it is also considerably sim- 
pler, since the incoming photon is a much better understood object than an incoming 
hadron. Aside from being a good testing ground of our understanding of perturbative 
QCD, it is also a probe of the structure of the target hadron. In fact, it has been 
often pointed out that heavy-quark photoproduction is a viable way to measure the 
gluon structure function in the proton 0. 

Radiative corrections to the single-inclusive photoproduction of heavy quarks have 
been first computed in ref. 0. The recent work of ref. || has confirmed the first 
computation, thus making the photoproduction cross section up to order 0(a em a 2 s ) 
a well-established result. From the next-to- leading-order computations the following 
facts have emerged. First of all, the photoproduction cross section receives more 
moderate next-to-leading corrections than the hadroproduction case. This result 
has improved the consistency of the data on charm production with the theoretical 
computation. In fact, before the radiative corrections were known, it was difficult 
to accommodate the experimentally observed hadroproduction and photoproduction 
cross sections with the same value of the charm quark mass, the first one requiring 
much smaller masses. 

A large amount of experimental information is available on photoproduction of 
heavy flavours^. Comparison between theory and experiments has not gone much 
further than the total cross section. This is a consequence of the fact that only charm 
production data have been available, and that the single- inclusive charm spectrum 
is strongly modified by non-perturbative effects. There is reasonable hope that, by 
looking at more exclusive distributions, we could learn more from photoproduction 
results. Modern fixed-target photoproduction experiments have the capability to 
study correlations between the heavy quark and antiquark. Furthermore, at the ep 
collider HERA, a large charm and bottom cross section is expected. It is clear, 
therefore, that in order to make progress in the physics of heavy-quark production, 
an exclusive next-to-leading-order calculation of the photoproduction cross section is 
needed. This may turn out to be useful both in charm production at fixed-target 
experiments and at HERA, and in bottom production at HERA. Since higher-order 
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corrections are moderate even in the charm case, it is possible that certain charm 
distributions may be used for QCD studies. 

In this paper we describe a next-to-leading-order computation of the doubly differ- 
ential cross section for the photoproduction of heavy-quark pairs. This computation 
follows closely the analogous work of refs. || and for hadroproduction of heavy 
quarks. Our result is implemented in the form of a "parton" event generator, which 
can be used to compute any distribution accurate to the next-to-leading order in the 
strong coupling constant. The problems arising from soft and collinear divergences 
are dealt with by generating appropriate sequences of correlated events, in such a way 
that the cancellation of collinear and soft singularities takes place for any well-defined 
physical distribution (i.e. distributions that are insensitive to soft and collinear emis- 
sion). The advantage of this method (developed for the first time in ref. 0) is that 
it does not require any artificial regularization of the cross section for producing the 
quark-antiquark pair plus a light parton. A detailed description of this method is 
given in ref. |5[]. In what follows we will describe the photoproduction calculation, 
with some emphasis on the differences with the hadroproduction case. 

This paper is organized as follows. In Section 2 we give a general description of 
the calculation. Some subtleties arise in the photoproduction calculation, which have 
to do with factorization scale choices. We discuss these problems in Section 3. Some 
phenomenological applications of our result have already been given in ref. |§, where 
a particular doubly differential cross section (of interest to the extraction of the gluon 
density from heavy-quark photoproduction data) is studied. In this work, we limit 
ourselves to the study of fixed-target photoproduction. More detailed studies of heavy 
quark production at HERA will be given in future works'! . In Section 4 we discuss 
the total cross section, and in Section 5 we discuss the differential distributions in 
fixed-target experiments. 

2. Description of the calculation 

The partonic subprocesses relevant for heavy-quark photoproduction at order 
a era a 2 s are the two-body process 



19 ^QQ 




-3- 



and the three-body processes 

19 -> QQg 

iq - QQq. (2.2) 
We will describe the two-body process in terms of the quantities 

s = ( Pl+P2 ) 2 

t = (pi - kx) 2 - m 2 = (p 2 - k 2 ) 2 - m 2 

u = (pi-k 2 ) 2 -m 2 = (p 2 -k 1 ) 2 -m 2 , (2.3) 

where p\ is the photon momentum, p 2 is the gluon momentum, and k\, k 2 are the 
momenta of the heavy quark and antiquark, respectively. We have p 2 = p 2 = and 
k\ = k\ = m 2 , where m is the mass of the heavy quarks, and s + 1 + u = (notice 
that the definition of t and u is not the conventional one). 

We will use dimensional regularization to deal with the divergences appearing in 
intermediate steps of the calculation. For this reason, we will need the expressions of 
phase spaces in d = 4 — 2e dimensions. The two-body phase space is given by 

2 2e /4ir\ £ 1 

* = F(T^(T) Ji-B^s^BM, (2.4) 

where (3 = y/1 — p, p = Am 2 / s and Q x is the angle between pi and k\ in the centre- 
of-mass system of the incoming partons. Therefore, 

* = -|(l-/3cos0i). (2.5) 
The three-body processes are characterized by five independent scalar quantities: 



s = 


(pi 


+ P2? 


tk = 


(pi 


-kf 


u k = 


(P2 


-kf 


11 = 


(Pi 


-hf 


Q2 = 


(P2 





(2.6) 
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where p\ is the photon momentum, p 2 is the momentum of the incoming parton, k± 
and k 2 are the momenta of the heavy quark and antiquark, respectively, and k is 
the momentum of the emitted light parton. We will often use the variable s 2 , the 
invariant mass of the heavy quark- antiquark pair, which is related to our independent 
invariants through 

s 2 = (k 1 + k 2 ) 2 = s + t k + u k . (2.7) 
It will be convenient to introduce variables x and y, where x = s 2 /s and y is the 

— * 

cosine of the angle between p~[ and k in the centre-of-mass system of the incoming 
partons. We have 

p < x < 1, -1 < y < 1 (2.8) 



and 



t k = ~(l-x)(l-y), u k = ~(l-x){l + y). (2.9) 



In the centre-of-mass frame of the QQ system, our four-momenta are given by 

Pl = pi (1,0, 0,1) 
p 2 = p\ (l,O,sin'0,cos'0) 
k = k° (1, 0, sin^', cos ■?/>') 



k\ = — — (1, f3 x sin 9 2 sin 6>i, f3 x cos 9 2 sin 9\, f3 x cos 9\) 



where 



k 2 = (1, — (3 X sin 6* 2 sin 9i, — (3 X cos9 2 sin^i, — f3 x cos#i), (2-10) 



_ s + t fc o s + u k _ i fe + M fc 
Pi ~~ o . — ) P2 ~~ : — ) ft — 



2^' ^ 2 2^' 2^ 



cos-0 = 1 sin-?/' > 

2P1P2 



P x = I/I-—. (2.H) 

V src 
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The two remaining independent invariants gi,?2 are given by 
9i = ~ S+ 2 k ( 1 ~ AkCOS0i) 

g 2 = fc (1 + /3a, cos 6*2 sinfli sin ^ + Ac cosffi cos ■?/>)■ (2-12) 

Now all invariants are expressed in terms of x, y, 9i, 9 2 and s through eqs. 
27TID and (111). 



The three-body phase space in terms of the variables x, y, 9i, 9 2 is given by 



,1-6 



where 



and 



d$ 3 = HNd^^—il -xf^Hl -y 2 )- e dysin- 2e 9 2 d92, (2.13) 

Z7T 



# = , F( ! / } r = 1 - -e 2 + 0(e 3 ) (2.14) 

r(l + e)r(l-2e) 3 V ; V ; 

N = ^r(l + c) (2.15) 



^ = r ,f £ v f— Y^/^sm-^xdcosMz. (2.16) 

1(1 — 6) \SX/ l07T 



Both 9 1 and # 2 range between and ix. 

We are now ready to compute the cross section for the real emission processes 
of eq. fl2.2|) . The technique is the same as that used in ref. ||. We begin with the 
subprocess 75 — > QQg. The cross section (in d space-time dimensions) is given by 

da$ = M^(s,t k ,u k , qi ,q 2 )d<S> 3 (2.17) 
M$( S , tM , qi) = ^-^^ ^^f, (2,3) 



where A$ is the invariant amplitude. The invariant cross section JA$?j has singular- 
ities in t k = and = 0, corresponding to soft (x = 1) and collinear (y = —1) gluon 
emission. No collinear emission from the photon line takes place at this order for the 
19 ~~ * QQg subprocess, and therefore M-jg is regular at y = 1. It can be shown that 
the leading soft singularity behaves like 1/(1 — a;) 2 , and that no double poles appear 
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in t k and u k . Therefore the function 

f jg (x,y,9 1 ,9 2 ) = 4t k u k M$ (s,t k , u k , q u q 2 ) (2.19) 

is regular for y = — 1 and x — 1 (the dependence of f ig upon s/m 2 is not explicitly 
shown). Using eqs. (|2.9|) we get 



A*y(M*,« fc ,gi,g 2 ) = g2 J 3 _ x)2(1 _ j/2 v ( 2 - 20 ) 



The three-body contribution to our cross section, including the phase space, is then 
given by 

da§ = HNd4 x)S ^dysm-^ 9 2 d9 2 (l ~ x)- l -^(l ~ yY^f^y^^^)- (2.21) 
We can now use the following expansions, valid for small e 

(!-*)—< = -£l,(l-,) + (_U 



2e V 1 — xjp \ 1 — x I . 

+C(e 2 ) (2.22) 

.(2a;)- 



l-y 2 )- 1 ' 6 = -[5(l + y) + 5(l-y)\- 



2e 



1 \ f 1 



H-O(c), (2.23) 



where the distributions in round brackets are defined according to the prescriptions 

V)(— ) = ftt 

p VI — x/p Jp 1 — x 

P \ 1 - x y . Jp 1 - x 

"•%)(rh) * - r 1 ^^^ (2.24, 

for any test function h(x). We define /3 = a/1 — p. The parameters p and u; should 
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be chosen within the ranges 



p < p < 1 , < w < 2. 



(2.25) 



The final results will not depend on the particular values chosen for p and u, but dif- 
ferent choices can lead to better convergence in the numerical programs, as discussed 
in ref. H. We obtain 



da$ = dai s ] + HNd® 



(x)< 



-l-e 



cfo/sin 2e 8 2 d8 2 



19 



X 



1 — X 



2e 



2tt 
logfl — X 



1 — X 



(l-l/ 2 )- 1 - 6 /^*, (2.26) 



where 



-l-e 



// A ^l 1 '— — cfy sin- 2e 6 2 d6 2 



2tt 



/9- 



-4e 



2e 



-6(1 -x) 



(l-y 2 )-^f jg (x,y,9 u 9 2 ). 



[2.27) 



The details of the calculation of the soft component of the cross section, dajfj, are 



given in Appendix A. Equation ( 2.27 ) can be explicitly integrated over x, y and 9 2 to 
obtain 



dai") = -HNd<$>, 



19 



'4ne 



- 1 - e p- 4e f\ s g \o 1 ) 



[2.28) 



where the function f!fj(9\) is given in Appendix A. 

We now expand (1 — y 2 )~ 1_e in the second term of eq. ( 2.26 ), observing that we 



only need the expansion up to order e°. As noticed above, M§j is regular at y — 1, 
and therefore the term proportional to 6(1 — y) gives no contribution. We get 



[2.29) 



where 



-l-e 



do-OH = Nd$F—— dysnT 2e 9 2 d9 2 
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2tt 



1 — X 



2e 



log(l — x) 
1 — X 



(2a;)" 



2e 



-6(1 +y) 



fi 9 (x,y,9x,9 2 ) 



(2.30) 
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and 



daifj = N-^^i3 x d cos 9 1 d9 2 dydx 



x 



1 



1 — X 



p L 



1 



1-3/, 



+ 



1 + 3/, 



f^{x,y,6 1 ,6 2 ). (2.31) 



The technique for the evaluation of the collinear limit of the invariant cross section 
is described in detail in ref. ||. The result in our case is 



* 1 e /2 \ ^(z) 



,y 4e Vu; 



1 — x 



2e 



log(l — x) 
1 — X 



/g->(Mi)> (2-32) 



where 

/(fWi) = 64^4^ (1 - x) 



1 — x 



1 — X 



X 



+ x(l — x) 



.Mg(xs )(?1 ). (2.33) 



Here 4 = a s /i 2e is the dimensionful coupling constant in d dimensions (the suffix 
(b) stands for bare), and M.}®(s, t) is the invariant cross section for ^yg — > QQ at the 
Born level, 



da® = M®(s,t)d$ 
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(2.34) 



The explicit expression for M.^ is given in Appendix A. The term in the square 
bracket in eq. ( [2.33j ) is, up to a factor 2C A , the gluon-gluon Altarelli-Parisi splitting 
function in d dimensions for x < 1, and the Born cross section is taken in d dimensions. 



With the usual definition 



- = --7e + log(47r), 
e e 



(2.35) 



we can rewrite eq. (|2.32|) as 



d< -> = -€1(1)' £m 



(6) 



1 — X 



X 



X + — + x(l - xY 



x 



2e / log(l - x) 
e I 1-x 



M { " ) Axs,q l )d^ ) . 



^2.36) 



We see that the 1/e divergence in the collinear term assumes the form dictated by 
the factorization theorem. According to this factorization theorem0, any partonic 
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cross section can be written as 



da ij (p 1 ,p 2 ) = Y^ / da k i(x 1 p 1 ,x 2 p 2 )T ki (x 1 )ri j (x 2 )dx 1 dx 2 , (2.37) 

kl J 



where 



Tij(x) = SijS(l - x) - \~Pij{x) + ^-Kij(x) + 0(a 

e In In 



(2.38) 



and da is free of singularities as e goes to zero. The collinear factors Tij(x) are 
usually reabsorbed into the hadronic structure functions, and only the quantities da k i 
will enter the physical cross section. The functions P%j{x) are the Altarelli-Parisi 
kernels. The functions K{j(x) in eq. (|2.38|) are completely arbitrary, different choices 
corresponding to different subtraction schemes. The choice Kij(x) = 0, to which we 
stick in the following, corresponds to the MS subtraction scheme^. 

Expanding eq. Q2.37Q order by order in perturbation theory, we find for our case 



1 Ctc 



da ig {piip 2 ) = da ig (p l ,p 2 ) + -—P gg (x)M\ ) g {xs : q 1 )d^ K 2 



(x) 



(2.39) 



where 



2C A 
2C A 



x 



X 



1 — X . 
H h x(l — x) 



+ 



X 

1 — X 



X 



+ x(l — x) 



+ 2nb 5(l - x) 

+ (2nb + 4C A \ogp)6(l-x) 



(2.40) 



and 



l\C A - AT F n u 

127T ' 



Here is the number of light flavours, and for N c = 3 we have 

1 

C 'a = 3, T F = — . 



(2.41) 



^2.42) 



The final expression for the short- distance cross section, after subtraction of the 
collinear divergences, eq. Q2.39|) , becomes 



da ig = da® + da^-) + d*$ + da<$ + da$ , 



(2.43) 
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where da$ and daifj are as given in eqs. (|2.34j) and ( |2.31|) , respectively, 



C A a< 



7T 



log ± + log | 



1 — X 



+ 2 



log(l — x) 
1 — X 



X 



X + 



x(l — x) 



A4jg(xs,<7i)d$£ 



(x) 
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da$ + -^-(27t6 + 4C A log^g( S ,i)^ 2 . 
7r e 



The scale /x, appearing explicitly in the expression of aig , has been set equal to a 
scale fi F characteristic of the subtraction of the singularity due to collinear emission 
from the incoming gluon, while a s is taken everywhere at the renormalization scale 
fi R (the problem of scale definitions and choices will be discussed extensively in the 
next section). The remaining singularities in daj^J are cancelled by the singularities 
in the virtual contribution to the cross section da^J . Therefore the quantity 



19 



(2.44) 
(2.45) 



da\f = do% + da% 



(2.46) 



is finite as e — > 0, and so is the full expression for da, 
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We now turn to the other three-body subprocess present at the a era a 2 s level, namely 
jq — > QQq. In this case, collinear emission takes place both from the photon and from 
the incoming light quark. On the other hand, there is no order a em a s contribution to 
heavy-quark pair production via jq fusion. Therefore, no soft singularity is expected. 
We explicitly checked that this is indeed the case. The three-body cross section is 
given by 



dat} 



i<i 



M$(s,t k ,u k ,qx,q 2 ) 



M i ^{s,t k ,u k ,q 1 ,q 2 )d^ 3 

1 1 V \A<r) 

2s2(l-e)N c A. I 79 

\ / ° spin, color 



(2.47) 
(2.48) 



where is the invariant amplitude. The invariant cross section A4^ has singu- 
larities for t k — 0, corresponding to collinear (y = 1) light-quark emission from the 
photon, or u k = 0, corresponding to collinear (y = —1) light-quark emission from the 
incoming quark. The function 



f iq (x,y,9 1 ,9 2 ) = 4t k u k M$(s,t k ,u k , q 1: q 2 ) 



(2.49) 
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is therefore regular for y = ±1, and vanishes as (1 — x) 2 for x — > 1. For this reason, 
the three-body cross section can be rewritten as 



da§ = HNd<S> { 2 x) ^—dysm- 2e 6 2 d9 2 



1 — x 



2e 



log(l — x) 
1 — X 



(l-y 2 )-^f, q (x } y } e 1 ,e 2 ).(2.50) 



Expanding (1 — y ) e in eq. ( |2.50|) we get 



(2.51) 



where 



da {c±) 
aa iq 



Nd<S> 



(*)■ 



-l-e 



2vr 



1 — x 



2e 



log(l — x) 

1 — X 



(2w) 



2e 



f 7g (x,y,9 1 ,9 2 ) 



(2.52) 



and 



(icr^ = N-^-^P^d cos 6id6 2 dydx 



1 — x 



1 \ / 1 

+ 



f iq {x,y,e u 6 2 ). (2.53) 



Performing the ?/ and 02 integrations in eq. (|2.52| ) we obtain 



dcr (c±) = -TV 

79 4e W 



d§\ 



1 — x 



2e 



log(l — x) 

1 — X 



4?^A)- (2-54) 



The explicit form of f^\x, 9i) can be obtained following ref. |J. The results are 

4 c g +) (x,^) = 32rta cm e 2 q s(l - x)P qj (x)M q b j (xs, q 2 ) , (2.55) 

where e q is the charge of the emitted light quark in electron charge units, M. q b q {s,t) 
is the lowest-order invariant cross section for qq — > QQ, 



d °m =M < S(s,t)d^ 2 , 



[2.56) 
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and 



P qi (x) 



x z + (l 



x 



2x(l — x)e 



(2.57) 



is the Altarelli-Parisi splitting function in 4 — 2e dimensions entering the probability 
of finding a quark in a photon, which is clearly equal to the splitting function of a 
quark in a gluon, up to a colour factor T F /N C , due to the A matrix in the gqq vertex. 
For f( c ~^ we obtain 



&-\x, 6x) = 327ra 5 s(l - x)P gq (x)M®{xs, Ql ), 



where 



1 + (1 - xY - ex 2 



x 



(2.58) 



(2.59) 



is the Altarelli-Parisi splitting function of a gluon into a quark in 4 — 2e dimensions. 
The subtraction of collinear singularities takes place as discussed for the previous 
case. We just give our final formulae, 



da iq = da^ q +) + da^ + dc$ , 

where da^ is given in eq. ( p.53|) , 

N c a em e 2 



(2.60) 



da {c+) 
aa iq 



2tt 



M q ^(xs,q 2 )d^ 



X 



2x(l -x) + (x 2 + (1 - x) 2 ) I log + log %- + 2 log(l 

V A* 7 2 



x) 



C F a s 
2tt 



X 



X + 



1 + (1 



x 



S LU 

log — + log - + 2 log(l - x) 



(2.61) 



(2.62) 



Notice that collinear emission from the photon is characterized by a scale /i 7 which 
is a priori different from the hadronic factorization scale jx F . The quantities dayj in 
eq. ( |2.43 ) and da^ in eq. ( |2.60| ) can be found in the literature^, and we did not 
need to explicitly evaluate them. The quantity da&J in eq. ( |2.43| ) was obtained from 
the authors of ref. 0. 

The analytical results presented here are implemented as a parton event generator, 
written in FORTRAN. The interested reader can obtain the code from the authors. 
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3. Scales in the photoproduction process 

Heavy-quark photoproduction differs from hadroproduction in the treatment of 
collinear singularities. In fact, when a light parton is collinearly emitted by the in- 
coming photon, the subtracted term is a signal from the non-perturbative region 
where the photon splits into quarks and gluons before interacting with the partons in 
the hadronic target. This fact is taken into account by inserting in the photon-hadron 
cross section a contribution in which the photon is formally treated as a hadron (the 
so-called hadronic or resolved photon component, to distinguish it from the point-like 
or pure-photon component, in which the photon directly couples with the partons of 
the hadronic target). The photon structure functions will also depend upon the mo- 
mentum scale /x 7 at which the collinear singularities of the photon leg are subtracted. 
Neither the point-like nor the hadronic components are separately independent of 
/x 7 , because the subtracted term in the point-like component is responsible for the 
redefinition of the photon structure functions in the hadronic component. 

Let us consider first the heavy-quark production process of an on-shell photon 
colliding with a hadron H at centre-of-mass energy yS. In order to clarify the role 
of the various scale dependences in the process, we write the 0(a em a 2 ) cross section 
in the following form 

Vq£\S) = J dx fj; H \x, fx F )a yi (xS,a s (fi R ), fi R , /j f , /j-y) 

i 

+ J dx ^ dx * /i 7) ( x i> A*7> ^F)fj H \ x 2, /i / F )o-i i (xix 2 S', a s (ji' B ), a4, (i' F ) 



'.J 



+0(a em a 3 s ) (3.1) 



with 



a. 



yt{s,a s (fj, R ),fj, R , fip,^) = a em a s (fi B )&-^ (s) + a em a 2 s (/j, R )a$ (s, /i R , /i F , /x 7 ) 
&ij(s,a 8 {n R ),n a ,H F ) = a 2 s {fi R )alf (s) + a 3 s {fi R )aff (s, fjt R , fi F , // 7 ). (3.2) 

Here fi R and fjf R are renormalization scales, /jl f and jjf F are factorization scales for 
collinear singularities arising from strong interactions, and factorization scale 

for collinear singularities arising from the electromagnetic vertex. If one wanted to 
extend eq. (|3.1|) to even higher orders, one should also include an explicit dependence 
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of the structure functions upon the renormalization scale. At the order we are con- 
sidering, the renormalization scale in the structure functions can be kept equal to the 
factorization scale, as is usually done. The left-hand side is independent of all the 
scales up to terms of order a em a|, provided the parton density functions obey the 
appropriate evolution equations. The hadronic and photonic parton densities obey 
the usual Altarelli-Parisi equations in /i F . In addition, the photonic parton densities 
have also an inhomogeneous evolution in // 7 , which, at the leading order, is given by 

d log /J,* 2tt 1 

where ei is the charge of the parton i in electron charge units. The compensation of 
the scale dependence takes place in the following way. The p, R scale dependence is 
compensated in the expressions for the partonic cross sections: the scale dependence 
of a s in the Born term a$ is compensated by the explicit scale dependence of the 
next-to- leading term o5 . A similar cancellation occurs in The dependence upon 
H F (fi' F ) cancels between the explicit dependence in the next-to-leading order term and 
the dependence in the structure functions convoluted with the Born terms. This holds 
independently for the two terms of eq. (|3.1j) . The dependence upon /i 7 cancels between 
the explicit dependence in the next-to-leading order component of the first term of 
eq. Q3.1]) and the /x 7 dependence of /j , as given in eq. (|3.3|) , multiplied by the Born 
level partonic cross section. In the commonly-used photon density parametrizations, 
/x 7 is usually kept equal to /i F , so that the term given in eq. ( |3.3|) becomes a correction 
to the usual Altarelli-Parisi equation (the so-called inhomogeneous term). Therefore, 
in our calculation, we use for consistency fi' F = /i 7 . In some cases, the inclusion of the 
hadronic component gives only a small effect, and will be neglected. In these cases 
we have chosen /i 7 = 1 GeV, which amounts to setting the photon structure function 
to zero at a scale of the order of a typical hadron mass. We have found that varying 
/i 7 between 0.1 and 5 GeV does not affect the results in a noticeable way. 

4. Total cross sections 

We begin our phenomenological study with the total cross sections for charm and 
bottom production. We will concentrate here on the analysis of the dependence of 
the total cross sections on the input parameters of the calculation, namely the choice 



x 



+ (l-x) 2 +0(a s ) 



(3.3) 
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of parton distribution functions (PDFs), the choice of quark mass and the choice of 
the values of the scales discussed in Section 3. A discussion of the ranges of masses 
and scales within which to explore the dependence of the cross sections can be found 
in ref. ||. 

The target will always be an isosinglet nucleon, N = (p + n)/2, and unless other- 
wise stated we will use the parton distribution set MRSD0& 

The default values of the charm and bottom mass will be 1.5 and 4.75 GeV 
respectively, and the default choices for fi F and fi R will be: 

\i R = m c , {i F = 2m c (4.1) 

for charm and 

fjL R = n F = m b (4.2) 

for bottom. The asymmetry in the default choice for the charm is related to the scale 
threshold below which PDFs extrapolations are not available, as explained in detail 
in ref. ||. As for /i 7 , we fix as a default // 7 = 1 GeV, as explained in Section 3. 

As an illustration of the reliability of the theoretical prediction we present in fig. [I] 
(fig. H) the leading and next-to-leading results for the total charm (bottom) cross 
section. The bands in figs. [I] and |2| are obtained by varying only the renormalization 
and factorization scales, everything else being kept fixed. 

Figures p] and |2| deserve some comments. First of all, the scale uncertainty associ- 
ated with the charm production cross section is significantly smaller at the next-to- 
leading order for beam energies above 200 GeV, while below 200 GeV the uncertainties 
at leading and next-to-leading order are similar. The residual uncertainty is much 
smaller than in the case of hadroproduction for comparable beam energies. The con- 
tribution of the hadronic component of the photon is always smaller than 5% of the 
total cross section for current energies. Uncertainties coming from the determination 
of the photon structure functions are therefore negligible with respect to others. 

The reduction of the variation band is even more pronounced in the case of bottom 
production, once next-to-leading order corrections are included. This is what we 
expect. For higher masses the value of a s is smaller, and the perturbative expansion 
becomes more reliable. Observe that the size of the leading-order band for the bottom 
cross section is not much smaller than the one for charm. This is due to the fact that, 
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as explained in ref. ||, we did not try to study the factorization scale dependence in 
the case of charm production. The reader should therefore remember that a further 
uncertainty should be added to the charm result, and that the band shown in fig. [I] 
is only an underestimate of the uncertainties involved in the computation of charm 
production cross sections. 

We now turn to all other sources of uncertainties, such as the structure function 
choice, the value of A QCD , and the mass of the heavy quarks. In tables [I] and |2| we 
give the cross sections for 7-nucleon collisions at various beam energies. The rates 
were obtained using a reference scale /io = for charm and /io = m& for bottom. We 
also show the effect of varying m c between 1.2 and 1.8 GeV, and m& between 4.5 and 
5 GeV. The scale fj, R was varied between /z /2 and 2/i . In the charm case, fi F was 
kept equal to 2/i , while for bottom we kept \i F = fi R . We verified that considering 
independent variations for the factorization and renormalization scale does not lead 
to a wider range in the bottom cross sections for the energies shown in the tables. 

The tables are broken into three blocks, each corresponding to a different choice 
of A QCD within the range allowed by the current uncertainties. The upper block 
represents the default choice relative to the MRSDO fit. The second and third blocks 



correspond to the sets discussed in ref. |16| for the nucleon. The values of A4 obtained 
in the fits of ref. |16| range from 135 to 235 MeV, corresponding to a range for A5 
between 84 and 155 MeV. This range for A5 is chosen because no good fit to deep 
inelastic data is possible outside that range in the context of ref. [[31J (i.e. with that 
choice of structure function parametrization, etc.). We have chosen instead the wider 
rangeB 100 < A 4 < 300 MeV, corresponding to 60 < A 5 < 204 MeV. Therefore, in 
order to take into account the full range of uncertainty associated with the value of 
Aqcd, we were forced to account only partially for the correlation beteween A QCD and 
the nucleon structure functions. 

For the charm cross section, as can be seen, the value of the charm quark mass is 
the major source of uncertainty. Differences between the extreme choices m c = 1.2 
and 1.8 GeV vary from a factor of 5 at 100 GeV to a factor of 3 at 400 GeV. Differences 
due to the scale choice are of the order of a factor of 2 at low energy, and at most 
50% at higher energy. A factor of 2 uncertainty also comes from the variation of A QCD 
within the chosen range. 

We also explored independently the effect of varying /i 7 . Differences are totally 
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negligible for values 0.1 GeV < /x 7 < 5 GeV, and are not included in the tables. 

All of these uncertainty factors are systematically a factor of 2 or more smaller 
than in the case of hadroproduction. Notice however one pathology encountered when 
combining the most extreme choices of mass (m c = 1.2 GeV), scale (fi R = m c /2) 
and A QCD (MRS235, A 5 = 204 MeV): the cross section in this case decreases in the 
region 100 GeV < E 1 < 600 GeV. This happens because for this particular choice of 
parameters the radiative corrections become negative in part of this energy range. 

Low-energy measurements of charm photoproduction cross sections favour a mass 
value of approximately 1.5 GeV. Previous comparisons with theory, however, were 
made using a fixed value fi R = 2m c . We expect that once the uncertainty on A QCD 
will be reduced, reliance on the next-to-leading-order calculation and the residual 
dependence on fi R should allow a determination of m c to within 100-200 MeV. 

The corresponding variations in the bottom case are smaller, in particular for 
energies sufficiently above the production threshold. At E 1 = 100 GeV, where we 
observe the largest uncertainty in the perturbative calculation, we should also expect 
large non-perturbative effects and therefore the perturbative prediction is not fully 
reliable. Notice also from fig. |2| that the contribution of the hadronic component of the 
photon represents, at low energy, a significant fraction of the total cross section. The 
reason is that, close to the threshold, the photon-gluon fusion process is suppressed 
by the small gluon density at large x, while production via the hadronic component 
of the photon can proceed through a light valence-quark annihilation channel. 

For completeness, we also give the contribution of the hadronic component of the 
photon in tables |3| and £|, evaluated using the photon structure function set ACFGP- 
mc of ref. |14| and the set LAC1 of ref. UlSf . As can be seen, at current fixed-target 
energies, this contribution is small. 

Our final prediction for the allowed range of charm and bottom production cross 
sections, including the full variation due to the scale choice, the value of A QCD and 
the nucleon structure functions is shown as a function of the beam energy in fig. [J[ 
for m c = 1.2, 1.5 and 1.8 GeV, and for m& = 4.5, 4.75, 5 GeV. The small contribution 
of the hadronic component is not included in the figure. 
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m c = 1.2 GeV 


m c = 1.5 GeV 


m c = 1.8 GeV 


E b /x R = 


m c /2 


m c 


2m c 


m c /2 


m c 


2m c 


m c /2 


m c 


2m c 


Nucleoli PDF set MRSDO, A 5 = 140 MeV 


60 GeV 


1.442 


0.937 


0.715 


0.533 


0.357 


0.275 


0.203 


0.138 


0.107 


100 GeV 


1.673 


1.205 


0.952 


0.746 


0.537 


0.427 


0.341 


0.248 


0.198 


200 GeV 


1.826 


1.506 


1.238 


0.979 


0.776 


0.640 


0.525 


0.414 


0.342 


400 GeV 


1.906 


1.753 


1.481 


1.160 


0.997 


0.844 


0.691 


0.582 


0.495 


Nucleoli PDF set MRS135, A 5 = 60 MeV 


60 GeV 


0.753 


0.627 


0.535 


0.312 


0.255 


0.216 


0.130 


0.105 


0.089 


100 GeV 


0.921 


0.807 


0.705 


0.443 


0.379 


0.329 


0.217 


0.183 


0.159 


200 GeV 


1.093 


1.016 


0.910 


0.602 


0.543 


0.483 


0.338 


0.300 


0.266 


400 GeV 


1.229 


1.190 


1.082 


0.740 


0.694 


0.629 


0.451 


0.417 


0.377 


Nucleon PDF set MRS235, A 5 = 204 MeV 


60 GeV 


2.794 


1.280 


0.886 


0.874 


0.467 


0.332 


0.304 


0.176 


0.127 


100 GeV 


3.012 


1.631 


1.179 


1.183 


0.700 


0.518 


0.502 


0.314 


0.236 


200 GeV 


2.923 


2.015 


1.530 


1.474 


1.002 


0.776 


0.751 


0.522 


0.409 


400 GeV 


2.714 


2.329 


1.830 


1.663 


1.280 


1.024 


0.956 


0.731 


0.592 



Table 1: Total charm cross sections (fib) in collisions, point-like photon contri- 
bution. Nucleon PDF set as indicated. 
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m b = 4.5 GeV 


m b = 4.75 GeV 


nib = 5 GeV 


E b li = 


m h /2 


m b 


2rrib 


m b /2 


m b 


2rrib 


m b /2 


m b 


2m,b 


Nucleoli PDF set MRSDO, A 5 = 140 MeV 


100 GeV 


0.081 


0.049 


0.031 


0.035 


0.020 


0.012 


0.013 


0.008 


0.004 


200 GeV 


1.092 


0.832 


0.615 


0.702 


0.525 


0.383 


0.445 


0.327 


0.235 


400 GeV 


3.969 


3.373 


2.760 


2.929 


2.469 


2.001 


2.163 


1.807 


1.452 


Nucleoli PDF set MRS135, A 5 = 60 MeV 


100 GeV 


0.083 


0.057 


0.038 


0.038 


0.025 


0.016 


0.016 


0.010 


0.006 


200 GeV 


0.962 


0.785 


0.618 


0.634 


0.509 


0.395 


0.413 


0.326 


0.250 


400 GeV 


3.279 


2.922 


2.511 


2.459 


2.174 


1.852 


1.845 


1.618 


1.366 


Nucleon PDF set MRS235, A 5 = 204 MeV 


100 GeV 


0.094 


0.055 


0.033 


0.039 


0.022 


0.013 


0.015 


0.008 


0.005 


200 GeV 


1.290 


0.951 


0.683 


0.823 


0.596 


0.422 


0.518 


0.368 


0.257 


400 GeV 


4.714 


3.909 


3.118 


3.459 


2.847 


2.251 


2.541 


2.074 


1.626 



Table 2: Total bottom cross sections (nb) in collisions, point-like photon contri- 
bution. Nucleon PDF set as indicated. 





m c = 1.2 GeV 


m c = 1.5 GeV 


m c = 1.8 GeV 


Eb /J-R = 


m c /2 


m c 


2m c 


m c /2 


m c 


2m c 


m c /2 


m c 


2m c 


Photon PDF set ACFGP-mc 


60 GeV 


0.050 


0.024 


0.015 


0.011 


0.007 


0.004 


0.004 


0.002 


0.002 


100 GeV 


0.096 


0.045 


0.027 


0.022 


0.013 


0.008 


0.007 


0.005 


0.003 


200 GeV 


0.203 


0.094 


0.056 


0.052 


0.029 


0.019 


0.017 


0.011 


0.007 


400 GeV 


0.376 


0.175 


0.105 


0.109 


0.059 


0.038 


0.039 


0.024 


0.016 


Photon PDF set LAC1 


60 GeV 


0.040 


0.022 


0.014 


0.010 


0.006 


0.004 


0.004 


0.002 


0.002 


100 GeV 


0.123 


0.053 


0.031 


0.019 


0.012 


0.008 


0.006 


0.004 


0.003 


200 GeV 


0.543 


0.205 


0.112 


0.075 


0.038 


0.023 


0.017 


0.011 


0.007 


400 GeV 


1.737 


0.647 


0.353 


0.292 


0.134 


0.080 


0.068 


0.037 


0.024 



Table 3: Contribution of the hadronic component of the photon to the total charm 
cross sections (/ib) in 7 A/" collisions. Nucleon PDF set is MRSDO, photon PDF set as 
indicated. 
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m b = 4.5 GeV 


m h = 4.75 GeV 


m b = 5 GeV 


E b n = 


m b /2 


m b 


2m b 


m b /2 


m b 


2m b 


m b /2 


m b 


2m b 


Photon PDF set ACFGP-mc 


100 GeV 


0.006 


0.006 


0.005 


0.003 


0.003 


0.003 


0.001 


0.001 


0.001 


200 GeV 


0.051 


0.059 


0.056 


0.035 


0.040 


0.038 


0.024 


0.027 


0.026 


400 GeV 


0.168 


0.189 


0.183 


0.123 


0.141 


0.138 


0.092 


0.107 


0.104 


Photon PDF set LAC1 


100 GeV 


0.005 


0.005 


0.004 


0.002 


0.002 


0.002 


0.001 


0.001 


0.001 


200 GeV 


0.054 


0.058 


0.053 


0.037 


0.039 


0.035 


0.024 


0.025 


0.023 


400 GeV 


0.163 


0.195 


0.187 


0.124 


0.146 


0.140 


0.095 


0.110 


0.105 



Table 4: Contribution of the hadronic component of the photon to the total bottom 
cross sections (nb) in 7 A" collisions. Nucleon PDF set is MRSD0, photon PDF set as 
indicated. 

5. Differential distributions 



In this section we present the results for the one- and two-particle differential 
distributions (the latter will also be referred to as correlations in the following) for 
charm and bottom production. Experiments can never obtain a monochromatic pho- 
ton beam, even though the energy of the photon can often be measured event by 
event. In order to retain a large statistics, therefore, differential distributions are 
most often presented integrating over the full photon beam spectrum, rather than at 
a fixed value of E 7 . While the total cross section is a slowly varying function of E 7 , 
features of the distributions can be affected by the convolution over the photon beam 
energy spectrum. A meaningful study therefore requires a similar convolution to be 
performed within the theoretical calculation. 

We have chosen to perform our phenomenological study using the photon beam 
energy spectra of the two experiments NA14/2 at CERN@ and E687 at FNAL@. 
The spectra were kindly provided to us by the two collaborations. Since the energy 
range of the two spectra are rather different, this choice can give an idea of the energy 
dependence of our results. The spectra are shown in fig. ^ and can be parametrized 
as follows: 



1 dN 

N dx 



C 



x ai (I - x) a2 + Rx bl (1 - x) b2 , (5.1 
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where x E^ / E max , E m i n <c E^ <C E max , 

E max = 440 GeV, E min = 125 GeV 
ai = 0.2294, a 2 = 2.9533 

i? = 7.5019 x 10 6 , 6i = 16.766, 62 = 11.190 (5.2) 

for E687, and 

E max = 400 GeV, E min = 20 GeV 
ai = 2.5987, a 2 = 4.7211 

H = -1.1862, 6i = 2.7133, 6 2 = 4.9516 (5.3) 

for NA14/2. The constant C is chosen to normalize the distributions to unity. The 
results we present are obtained by convoluting the theoretical distributions with the 
above beam shapes. We will refer to these, respectively, as NA14 beam and E687 
beam. 

Needless to say, only a detailed simulation of the detector acceptances and efficien- 
cies can allow a complete comparison of our results with the experimental findings. 
Therefore one should take the results presented here as indicative of the most relevant 
features of the next-to-leading-order calculation. Additional non-perturbative effects 
such as fragmentation and intrinsic momentum of the partons inside the hadrons will 
also be discussed at the end. 

As in the case of the total cross sections, the distributions will be calculated using 
the parton distribution set MRSDOU for the nucleon, unless otherwise stated. The 
default values of the masses will be m c = 1.5 GeV and m& = 4.75 GeV. The default 
values of the factorization and renormalization scales \i F and \x R will be: 



for charm, and 



for bottom, where 



— Ho, A*ii — A*o 



Hq = Jm 2 +p 2 



(5.4) 

(5.5) 
(5.6) 
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for the one-particle distributions and 

Mo = \/^ + ^±^ (5.7) 

for correlations. Here p T and p T are the transverse momenta of the heavy quark and 
antiquark, respectively. In the case of correlations, a further ambiguity arises in the 
choice of the scale (described in detail in ref. ||), which has to do with the freedom 
of choosing the same or different values of the scales for the three-body event and for 
the corresponding counter-events. The results presented here will always follow the 
approach of recomputing the scale for the counter-events. 



5.1. Charm production 

We begin with one-particle inclusive rates. On the left side of fig. [| we show the 
inclusive p T distribution of c quarks in the case of the E687 and NA14 photon beams. 
The solid lines represent the full next-to- leading-order result. The dots give the 
leading-order contribution rescaled by a constant factor. A slight stiffening of the p T 
distribution after radiative corrections is observed. 

On the right-hand side of fig. [| we show the inclusive x F distribution for the full 
next-to-leading-order calculation, superimposed onto the rescaled Born result. Our 
definition of x F is 

x F = -L (5.8) 

where CM refers to the centre of mass of the target and the tagged photon beam. In 
other words, this centre of mass has a boost with respect to the laboratory frame, 
which depends on the energy of the photon responsible for the interaction; py is the 
momentum projection on the beam direction in the centre-of-mass frame, and -Ecm is 
the centre-of-mass total energy. The fraction of c quarks produced with positive x F is 
larger than 90%, due to the hardness of the photon probe. Notice however that next- 
to-leading-order corrections induce a softening of the distribution, due to processes 
where a photon splits into a light-quark pair and interacts with a light quark from 
the nucleon. We verified that the inclusive p T spectrum does not change in shape 
if we restrict ourselves to quarks produced at positive x F , which is the region where 
data are usually collected by the experiments. 
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Next we consider cc correlations, starting from the invariant mass of the heavy 
quark pair, which is shown in fig. ^. As before, the continuous lines represent the full 
result of the next-to-leading-order calculation, while the dots are the rescaled Born 
results. The lower lines are obtained by imposing x F > for both quarks. Notice 
that while at the inclusive level most of the c's have positive x F , a large fraction of 
pairs with large invariant mass have either the c or the c at x F < 0, as one should 
expect. 

The difference in rapidity between the quark and the antiquark and the x F of 
the pair are shown in figs. [7] and |8|. A slight broadening at next-to-leading order is 
observed for the rapidity correlation, while a dramatic change is observed at next- 
to-leading order in the case of the pair x F distribution. This dramatic change can 
be traced back to the particular kinematics of heavy-quark photoproduction. In the 
Born approximation, the x F of the heavy-quark pair is simply related to its invariant 
mass 

M 2 - 

X QQ = 1 __9Q (5 . 9) 

and is therefore peaked at xj® near 1. At the next-to-leading level, it is the x F of the 
QQg system that will be peaked near 1. This means that the x F of the heavy-quark 
pair will be markedly degraded. The contribution of fusion is even softer, because 
of the contribution of the photon splitting into a light quark-antiquark pair. 

We now consider correlations that are trivial at the leading order, namely the 
difference in azimuth and total transverse momentum of the quark pair, p®Q. At 
leading order these distributions are delta functions centred respectively at A0 = ir 
and p®Q = 0. Higher-order real corrections such as gluon radiation or gluon splitting 
processes smear them. We plot these distributions in fig. ||. Even after the inclusion 
of higher-order effects, the azimuthal correlation shows a strong peak at A0 rs tt. 
Likewise, the p®Q distribution is dominated by configurations with pQQ smaller than 
1 GeV. In both cases, the tails are higher at the higher energy of the E687 experiment. 

For all the previous distributions, changing the values of the charm mass and 
renormalization scale fi R results in large differences in rates but small and easily 
predictable shape modifications. The pattern of these changes is similar to what is 
observed in the case of fixed-target hadroproduction® . 
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5.2. Bottom production 

The differential distributions for bottom production are shown in figs. |H] to 0. The 
cross section at NA14 is very small. We nevertheless include their distributions for 
completeness. 

The p\ (single-inclusive) and the Mqq distributions are broader than the cor- 
responding distributions for charm production. They are, however, narrower than 
would be expected on the basis of simple scaling arguments. This is because b pro- 
duction at fixed-target energies is still too close to the threshold, and thus constrained 
by phase-space effects (this can also be noticed from the strong energy dependence of 
the shape of the curves). For the same reason, the Ay distribution is narrower than 
in the charm case. 

The p®Q and the A0 distributions are much narrower, as a consequence of the 
smaller value of a s at the b mass, and of the previously mentioned phase-space con- 
straints. The x F distribution is softer for bottom than for charm, and in particular 
the fraction of b quarks with negative x F is larger. This is because a harder parton 
from the nucleon is required to reach the energy threshold for the creation of a b pair. 



5.3. Higher-Order Corrections and Hadronization 

The results discussed so far were obtained with a purely perturbative calculation 
limited to the next-to-leading order. In the case of charm quarks, the dependence 
of the cross section on the renormalization and factorization scales indicates that 
higher-order corrections might be large. Nevertheless, the stability of the shapes 
under inclusion of the next-to-leading-order terms suggests that no significant changes 
should be expected in the differential distributions when yet higher-order terms are 
included. This is not necessarily true of possible non-perturbative corrections, such 
as the intrinsic p T of the initial-state partons, hadronization and fragmentation. In 
particular, in the regions of phase space close to threshold, or for p® Q — > 0, we should 
expect significant corrections. 

In our previous study of heavy-quark correlations in hadropro duct ion @ we ex- 
plored these effects using different phenomenological models. In particular, we con- 
sidered the parton shower Monte Carlo HERWIG@ to simulate both the backward 
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evolution of the initial state and the formation of charmed hadrons. We found that the 
backward evolution of initial-state gluons gave rise to an artificially large intrinsic p T 
of the gluons themselves ((&t) = 1.7 GeV). This large intrinsic p T would significantly 
broaden the A<p correlation, in addition to softening the inclusive p T distribution. 
We also found that, as expected, the colour singlet cluster formation and decay into 
hadrons leads to a large colour drag in the direction of the hadron beam. The com- 
bination of the intrinsic p T smearing, of the colour drag, and of the decay of unstable 
charmed hadrons, resulted in inclusive p T and x F distributions of D mesons that are 
slightly harder than those resulting from the purely perturbative calculation. The 
A(j) correlations remain significantly broader than described by next-to-leading-order 
QCD. Both these results are supported by current experimental evidenced. 

In ref. we were also able to parametrize the effect of the intrinsic p T by using 
a Gaussian smearing, while we argued that there is no solid basis for the application 
of the Peterson formalism© to describe the fragmentation in the large-x F region. In 
this section, we repeat the analysis in the case of heavy-quark photoproduction, for 
the case of charm production with the E687 photon spectrum. 

Figure [13| shows a comparison of inclusive p T and x F distributions obtained from 
the perturbative calculation (solid lines), from HERWIG before hadronization (dot- 
ted lines), and from HERWIG after hadronization (dashed lines). In the case of the 
inclusive p T distribution we find that, as observed in hadropro duct ion, the effects of 
intrinsic p T and of hadronization tend to respectively harden and soften the distribu- 
tion. The net result is a softening of the shape. This is due to the fact that only one 
initial-state parton (namely the gluon) can acquire a transverse p T in photoproduc- 
tion, contrary to hadroproduction where we have two gluons in the initial state and 
the effect is enhanced. 

The most dramatic effects are however observed in the inclusive x F distribution. 
Hadronization effects heavily suppress the production of charmed hadrons at large 
x F , which was favoured at the purely perturbative level. This is perhaps surprising, 
because it would be expected that for heavy quarks produced in the photon frag- 
mentation region we should be able to obtain the correct distribution by using the 
perturbative calculation, convoluted with the fragmentation function for the c-quark 
fragmenting into a D meson. This is not necessarily true. As in the hadroproduction 
case, there are no theoretical reasons to support this possibility. The fragmentation 
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function formalism is in fact applicable only when the elementary production pro- 
cess takes place at energies much larger than the mass of the heavy quark. This 
is the case, for example, in the production of heavy flavoured mesons in e + e~ col- 
lisions, or in hadroproduction and photoproduction at high transverse momentum. 
The factorization theorem states that the same fragmentation function, evolved to 
the appropriate scale, and convoluted with the perturbative calculation of the par- 
tonic subprocess, should describe all these processes. The x P spectrum, instead, is 
not really characterized by a high-energy elementary process. The heavy-quark pair 
is produced with a relatively small invariant mass, and the large x F region is reached 
when the production angle in the heavy-quark centre of mass is small. Under these 
circumstances, non-perturbative effects (other than the fragmentation effects) could 
also take place. For example, the heavy quark could feel the dragging of the heavy 
antiquark, or of the beam remnants. 

Figure |T3| shows the HERWIG results for the charm pair p®® and A<f) distribu- 
tions. While the A0 correlation is still significantly broader than that calculated at 
the perturbative level, the distribution is more peaked than the one evaluated from 
HERWIG in the hadroproduction case. This results from the smaller effect of intrinsic 
p T present in the photoproduction, as already mentioned above. Existing data [BE! 
agree qualitatively with this result. 

In fig. [15|, finally, we show the charm pair p®® and A</> distributions obtained 
by giving a random intrinsic transverse momentum to the incoming gluon, with a 
Gaussian distribution!! , for different values of As can be seen, the choice (h?) = 
3 GeV 2 reproduces quite closely the HERWIG result for the A(f> correlation. As 
already discussed in ref . , we expect such a large intrinsic p T to be a pure artefact 
of the Monte Carlo. If the data were to confirm the existence of broad A0 correlations 
as shown in fig. [TJ], it might be interesting to try to justify theoretically on a more 
solid basis the possibile existence of such a large intrinsic p T for the gluons inside the 
hadron. 

5.4. Effect of the photon hadronic component 

As already seen in Section 4, the effect of the hadronic component on the total cross 
section, predicted using standard photon structure functions, is generally small. The 
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question now arises whether its effect on the shape of distributions is also negligible. 
We have studied this problem using the next-to-leading calculation of ref . || , together 



with the parametrization of the photon structure function set ACFGP-mc of ref. [14 



and the set LAC1 of ref. |{L5 |, at a fixed photon energy of 230 GeV. We have found 
that the only important modifications occur in the x F distribution for a single heavy 
quark, and for the x F of the pair. In the inclusive x F distribution for charm, we 
find that the hadronic component becomes comparable to the point-like term for 
x F pa —0.3, and for smaller x F it remains of the same order. For b production, we find 
instead that the hadronic component becomes comparable to the point-like term for 
x F pa —0.7, a region in which the cross section is several orders of magnitude below 
the peak value. 

For the x F of the pair, the effect of the hadronic component is more pronounced. 
This is due to the fact that the point-like component is concentrated near x F Q = 1. 
On the other hand the hadronic component is distributed in the central region (we 
find that for charm it peaks at x F pa —0.15 with the LAC1 set, and at x F ~ with the 
ACFGP-mc set). Its contribution in this region is therefore of the magnitude given 
in tables [| and Observe, however, that (as discussed in the previous subsection) 
hadronization effects do spread out the x F distribution of the pair. It is unlikely, 
therefore, that one can use these distributions to make statements about the hadronic 
component of the photon. 

We conclude, therefore, that for all practical purposes, the hadronic component 
can be neglected altogether in the fixed-target experimental configurations of present 
interest. 



6. Conclusions 

We have performed a calculation of next-to-leading-order QCD corrections to 
heavy-quark production in photon-hadron collisions. Our calculation improves over 
previous results, in that it can be used to compute any distribution in the heavy 
quark and antiquark, and possibly in the extra jet variables. 

We have presented a phenomenological study of total cross sections, single-inclusive 
distributions and correlations, for charm and bottom production at fixed-target en- 
ergies. 
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A detailed study of all the theoretical uncertainties of the calculation, namely 
those due to an independent variation of heavy-quark mass, factorization and renor- 
malization scales and A QCD , has been performed for the total cross sections. We found 
that the next-to-leading contribution is less important here than in the hadroproduc- 
tion case. Therefore, we expect the full result to be more reliable for photoproduction 
than for hadroproduction of heavy flavours. 

For single and double differential distributions, in the case of charm production, we 
always find that non-perturbative effects could be important. For the p T distribution, 
a possible description of the non-perturbative effects is given via the introduction of an 
intrinsic transverse momentum for the incoming gluon, which tends to stiffen the p T 
distribution, together with a fragmentation function similar to the ones used in e + e~ 
physics. This seems to give a reasonable description of the p T and A0 distributions 
in both photoproduction and hadroproduction. 

In the case of the x F distribution, it is very difficult to describe the non-perturbative 
effects in a simple way. From Monte Carlo studies, we conclude that it is likely that 
colour-dragging effects prevail, thus making it difficult to give a homogeneous descrip- 
tion of the hadroproduction and photoproduction x F distribution. 
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Soft limit of the real amplitude 

In this appendix we present the calculation of the amplitude for the process 7(7 — > 
QQg in the limit when the momentum of the emitted gluon tends to zero. This is 
the only case of interest, because the analogous limit for the process 75 — > QQq gives 
a trivial result. Momenta and colour and polarization indices are assigned as follows: 

7(Pi,/-0 + 9(p2,v,a) -> Q(h,i) + Q(h,j) + g(k,p,b). (.1) 

The main difference with respect to the heavy-quark hadroproduction case is that 
the gluon cannot directly couple to the incoming photon. In principle we are then 
left with three potentially singular diagrams when k — > 0, namely the diagrams in 
which the outgoing gluon is emitted by the incoming parton or by the heavy quark 
or antiquark. We indicate the amplitude for the process 7^ — ► QQ with 

where the momenta and indices are as in eq. ([!]). It then follows that the contribution 
of the diagram in which the gluon is emitted from the outgoing heavy quark is given 
by 

9s>hi x k + ki + m 

When the gluon is emitted from the antiquark we have instead 

9s>Hj- n w a k + h + m 

— u{hl) (k + k 2 y - 7 ' vih) ■ G4) 

Finally, when the gluon is emitted from the incoming gluon leg the contribution is 

crA 

^ fa _ k )2 ~ k > k ~ 6 (*0 A^ ;ij v(k 2 ) , (.5) 

where 

G t%» 3 (vii qa) = f abc tw»(9i - + sw 3 (?2 - g 3 ) w + Wife - (- 6 ) 



is the quantity that appears in the QCD three-gluon vertex. We can now evaluate 
the k — > limit in eqs. ([3|), (Q) and ([5]). The soft limit of the three-body real 
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amplitude is 



A1 ( k - °) 



2 fc- jfex 



+Wsf 



k-p 2 



(•7) 



At the lowest order, the amplitude for the two-body process jg — > QQ can be written 



as 



Using the identity 



eq. ([7|) becomes 

(* -> o) = 



fi(*i)A^«(A; a ) = A^B^. 



(.9) 



(A 6 A a ) 



ij \k ■ k\ k ■ p 2 



( A " A '), f 



ij \k ■ k 2 k ■ p 2 



(.10) 



Squaring and summing over initial and final degrees of freedom we obtain 



A% (k - 0) 



9$ 



Tr (A 6 A a A a A fe ) ((hh) - 2{p 2 k 1 ) - 2(p 2 k 2 ) + (k 2 k 2 



2Tr (A A" A A") (hk 2 ) - (Mi) - (p 2 k 2 ) 



E (B^) f , 
(.11) 



spin 



(.12) 



where we introduced the eikonal factors defined by 

/ j = v ■ w 

v ■ k w ■ k ' 

and we used the masslessness of the incoming parton, p\ = 0. Evaluating the traces 
Tr (A 6 A a A a A fc ) = 16T F D A C F , Tr (A 6 A a A b A a ) = 16T F D A (c f - -C^ , (.13) 
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inserting the proper flux factor, and averaging, we finally have the squared three-body 
amplitude in the soft limit 



M {s) 



gl 



C F \{k x kx) + {k 2 k 2 



2[C F - l -C A )(kyk 1} 



-C A {p 2 h) + (p 2 k 2 ) 



M {b) . 



;-i4) 



For QCD with three colours T F = 1/2, D A = 8, C F = 4/3 and C A = 3. Here 

1 1 



73 2s4D A 



Tr (\ a X a ) (B^) 1 " 



.15) 



spin 



is the Born amplitude squared for the two-body process jg — > QQ. By direct calcu- 
lation we obtain 



19 



e 2 n 2 T 
sut 



(f + u 2 -es) (l-e)+4m 2 s (1 



m 2 s 



ut 



;.i6) 



where e Q is the heavy-quark charge. Equation (|. 14|) can now be integrated over three- 
body phase space. Using the notations of Section 2 we have 
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dxdyd6 2 5(l - x) (l - y 2 ) 1 'sin 2e 9 2 [At k u k M 
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dxdyd6 2 5(l -x)(l- y 2 ) * 'sin^^ [U k u k M^] . (.17) 



The only non-trivial point of the calculation is the integration of the eikonal factors 



I{vw) = J dyd6 2 (l-y 2 ) 1 e sin- 2 ^ 2 [4t 



k u k ivw 



\x=l 



'It 



The expression for f( s J (O^ is obtained by formal substitution of each (yw) eikonal 
factor in eq. ( |.14[ ) for the corresponding integral I( vw ). We give here all the eikonal 
integrals, thus correcting some misprints of the analogous formulae in ref. ||. 



'(PlP2) 



2 W 

is— 

e 



ASTT 



-2eLi 2 1 + 



-t 
m 2 

21 



~t 6 , 2 

lo S — + lo S — ~ g g 



8{l-0). 



2eLi 2 1 + 



1 + 
1-0 

2t 



8(1+0), 
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^(P2fel) 
-^(pife) 

hk2k 2 ) 
hkik 2 ) 



-2elog 



-2t 



log- 



-2t 



8(1 + 13) *s(l-(3) 

J(pi*l)(* U ) 

I(p2k\) 

I{p\ fei ) 

- 1 + X±f) 



= /, 



(fcifei) 



8s 1 - 



p\ 7T 



27/5 



\ 1 + P ( 2(3 
log + e I Li 2 - r - Li, 



1-/9 



+ 



-2/j N 



.19) 
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Figure Captions 

Fig. 1: Total cross section at leading and next-to- leading order for charm produc- 
tion in 7iV collisions as a function of the beam energy. We plot the range 
of variation for the scale changes indicated in the figure, for m c = 1.5 GeV. 
MRSDO parton distribution set. The dotted line represents the contribu- 
tion from the hadronic component of the photon, evaluated using the set 
ACFGP-mc. 

Fig. 2: Total cross section at leading and next-to- leading order for bottom produc- 
tion in 7iV collisions as a function of the beam energy. We plot the range of 
variation for the scale changes indicated in the figure, for = 4.75 GeV. 
MRSDO parton distribution set. 

Fig. 3: Total cross section for charm and bottom production in 7iV collisions as 
a function of the beam energy. We plot the overall range of variation for 
changes in the parameters as discussed in the text, each band representing 
the result for a specified value of the quark mass. 

Fig. 4: Photon beam energy spectra for NA14/2 and for E687. 

Fig. 5: Charm inclusive p T and x F distributions in 7iV collisions with the E687 and 
NA14 photon beam energy spectrum. 

Fig. 6: Invariant-mass distribution of charm pairs produced in 7^ collisions with 
the E687 and NA14 photon beam energy spectrum. 

Fig. 7: Rapidity correlation and x F distribution for charm pairs produced in 7^ 
collisions with the NA14 photon beam energy spectrum. The lower curves 
are obtained requiring both quarks to have x F > 0. 

Fig. 8: Same as fig. 0, but for the E687 photon beam spectrum. 
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Fig. 9: Charm pair p™ and azimuthal correlation in <yN collisions with the E687 
and NA14 photon beam spectrum. 

Fig. 10: Bottom inclusive p T and x F distributions in >yN collisions with the E687 
and NA14 photon beam spectrum. 

Fig. 11: pQ® and invariant-mass distributions of bottom pairs produced in ^yN col- 
lisions with the E687 and NA14 photon beam spectrum. 

Fig. 12: Bottom pair correlations in 7^ collisions with the E687 (solid lines) and 
NA14 (dashed lines) photon beam energy spectrum: azimuthal correlations 
(left side), x F of the pair (right side) and bb radipity difference (left inset). 

Fig. 13: Comparison between HERWIG and the 0(a em al) result (solid) for inclusive 
distributions of charm with the E687 photon beam. For HERWIG we plot 
the variables relative to the charm quark before hadronization (dotted line) 
and relative to stable charm hadrons (dashed line). 

Fig. 14: Comparison between HERWIG and the 0(a em a 2 s ) result (solid) for inclusive 
pQ Q and A0 distributions of charm with the E687 photon beam. Different 
line patterns are explained in the previous figure's caption. 

Fig. 15: Effect of a non-perturbative px kick for the incoming parton in the nu- 
cleon, compared with the 0(a em a 2 s ) effect. The curves with a px kick are 
obtained with the Born cross section with MRSDO structure functions, sup- 
plemented by a random px kick on the incoming parton. The NLO curves 
are obtained with the same structure functions and were rescaled to the 
same normalization as the other curves. 



